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1 Introduction 

Group behavior has received much attention as a test case of self-organization. 
There has been much written in recent years to investigate interactions within 
groups of agents. These agents can be animals moving in an interactive way, such 
as birds, but can also refer to situations such as people driving in traffic. The 
models that describe these interactions are able to reproduce different structures 
and patterns relating to the movement and interaction of the agents involved. 
The advantages and necessities of this type of analysis in any complex biolog- 
ical, technological, economic, or social system are important and far-reaching. 
Each model that we will discuss describes interaction between agents. 
Focusing on animal groups, such systems received attention not only from biol- 
ogists and ethologists, but also from physicists, mathematicians, and engineers. 
This interest has produced a huge amount of literature, which is well docu- 
mented and reviewed: see, for instance, [ZlEllTl] and [3]. Loosely speaking, the 
basic idea behind these works is that complex collective behavior arises from 
simple interactions among close animals. 

This paper considers a simple case of five birds in V formation, flying at 
constant velocity. The aim is to study the feasibility of this solution and its 
stability for three models proposed in the literature. 

The first model we will analyze was proposed in [J. It is an agent-based model, in 
the sense that each agent is singularly considered. The model is also leaderless, 
meaning that all agents behave following the same set of rules and their behavior 
is not imposed by others. The model is purely based on attraction-repulsion 
infractions between group mates. Attraction allows the group to be formed 
and stay tight, while repulsion allows for collision avoidance between group 
mates and keeps them well spaced. This model does not combine attraction 
and repulsion with velocity alignment, unlike many other models that have 
been proposed. The model also assumes that each agent interacts with a limited 
number of group mates that are in a suitable pre-defined sensitivity zone. The 
authors were able to produce cluster formations, line formations and V-like 
formations by altering the attraction and repulsion parameters of their model. 
The second model was proposed by Cucker and Smale in [6] and additional 
results are available in [10]. The model is of second order, in the sense that 
the dynamics is given for position and velocity of each animal. More precisely. 
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each animal adjusts its velocity to the mating ones, which are in its "interaction 
network". The adjustment is done towards an average value of the various 
different velocities. Therefore this approach is quite related to the theory of 
consensus. 

The third model is again a consensus-type model proposed by Olfati-Saber in 
[5]. In this case, beside the adjustment of velocity towards an average of mating 
birds, there is a term which forces animal to keep a desired distance d. The 
latter is thought as a suitable distance which taken into account both attraction 
among birds and repulsion for avoidance of conflicting trajectories. 

For each model, we first show that the V formation is an admissible solution 
(for the first model we need to accurately choose parameters). Then we study 
stability of the linearization along the given V formation trajectory. The result 
is that the linearization of the first model gives rise to negative or zero eigenval- 
ues, while the other models presents also positive eigenvalues. We can conclude 
that the second and third model are exponentially unstable. 
In other words, despite the fact that these models can represent V like forma- 
tions, they are not able to represent stable V like formations. This indicates the 
need of further modelling efforts. 

2 Animal groups models 

Let us start with the model proposed in 0] , which we refer to as the CFP model. 
The animals in the model are represented by point particles, which have simple 
continuous-time dynamics. Given N G N, for alH = 1, . . . , and t £ M>o, let 
Xi {t) G represent the position of the «-th animal, whose evolution is described 
by the differential equation 



where x(t) is the vector {xi{t), . . . ,XNit)). 

The velocity Vi{x) is the sum of two contributions, expressing the effects of 
attraction and repulsion. 



In more detail, each of these contributions depends on the relative position of 
the other animals. 



Xi{t) = Vi{x{t)), 



(1) 



v^{x) = vt{x) 






The following definitions have been used. 
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• The function /a : R>o — >■ K>o (resp., /r : ]R>o — ^> R>o) describes how each 
animal is attracted (resp., repelled) by a neighbor at a given distance, 
assuming || • || denotes the Euclidean norm in R^. 

• The attraction neighborhood A" (resp., the repulsion neighborhood TZ") 
is the set of the n animals closest to the i-th one, which are inside the 
attraction (resp., repulsion) sensitivity zone. 

The above model is very general, and we need to specialize it by choosing 
the interaction functions and the shape of the sensitivity zones. We make the 
following assumptions. 

Al. The functions /a and /r are assumed to be 

M\\Xj - X.^\\) = F^\\Xj - X.^\\ , fr{\\Xj - X^\\) = -, 

\\Xj Xt\\ 

where -Fa and F^ arc two positive constants. 

A2. The sensitivity zones are depicted in Figure [T] and illustrated as follows. 
Let the center point be the animal's position, and let the horizontal axis 
(arrow-headed) represent the direction of motion. Attraction is active in 
a frontal cone whose width is given by the angle G (0, 360°] (dashed 
line). Repulsion is active both inside a disk of radius i?sr > {short-range 
repulsion) and in a frontal cone of width S (0°, 360°] (solid line). We 
stress that the sensitivity zones do not necessarily coincide with the visual 
field of the animal. They rather represent the zones which attraction and 
repulsion are focused on. 

A3. The speed of each animal \\vi\\ is bounded from above by a constant Wmax- 

The above assumptions result in the system 

Mt) = Fa 5: {X, - -f^y: ti'i - (2) 

Some remarks are in order. 

Rl. The definition of the interaction neighborhoods Af and 7?." allows to have 
an a priori bound on the number of effective neighbors, and therefore 
on the sensing and "computational" effort which is required for each an- 
imal. This fact, which copes with animals' intrinsic limitations, has been 
experimentally observed in biology, for fish [l] and birds [2]. The latter 
paper calls this neighborhood definition topological, as opposed to metric 
definitions, based on distance only. 

R2. Our assumption of unbounded sensitivity regions does not intend to imply 
that animals sensing capabilities extend on an unlimited range, but rather 
that group dynamics happen in a relatively small area. 
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Figure 1: The shape of the sensitivity zones. 



R3. If = Q!r = 360° the parameter Rg^ has no effect, and the interaction is 
completely isotropic as in the simulations presented in [2j[5]. Note that, 
even if there is no preference for any specific direction, the limitation of 
the number of considered neighbors makes the interactions not reciprocal, 
i.e. the fact that the i-th animal interacts with the j-th does not imply 
that the j-th interacts with the i-th. 

R4. By specializing the functions /a and /r, one can obtain various interaction 
models. For instance, it is natural to assume that /r be decreasing and 
lims^.oo /r(s) = 0. Similar considerations are valid also for some features 
introduced in other models, such as a neutral zone around animals or 
a hierarchical decision tree which allows the repulsion force to have the 
priority over the attraction force. 

The second model we will analyze is the Cucker-Smale model. It is inspired 
by the Cucker-Smale flocking model, and attempts to analyze animal behavior 
under hierarchical leadership. The goal of this flocking study is to be able 
to interpret, model, analyze, predict, and simulate various flocking or multi- 
agent behavior. The Cucker-Smale flocking model includes an interaction weight 
function which will allow us to analyze the birds in our formation by defining 
the weight function appropriately. In a model of hierarchical leadership, the 
agents are ranked in some way and are therefore not independently motivated 
in their movement - that is, lower-rank agents are led by agents of higher ranks. 
Using these assumptions, the paper goes on to discuss various flocks using these 
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assumptions. 

The Cucker-Smalc flocking model which is given by the autonomous dynamical 
system: 

1 = EjeL(j) ajj(a;)(vj - Vi), i = 1, . . . i > 

This model differs from the previous model in that the motion in each bird is 
defined in terms of Xi{t) and Vi{t). In other words, this is a second order model, 
assuming acceleration driven by "social forces". In this case L{i) denotes the 
subgroup of animals that directly influence agent i. 

The third model we will analyze is from [5] . It presents a theoretical frame- 
work for design and analysis of distributed flocking algorithms. In this paper, 
three flocking algorithms are presented: two are for free-flocking and one for 
constrained flocking. As in the Cucker-Smale model, this model compensates 
for velocity alignment. Their first algorithm generally leads to fragmentation, 
whereas the second and third algorithms both lead to fiocking behavior. How- 
ever, unlike the second model, this paper asserts that flocks do not need leaders. 
The paper also includes various simulation results that demonstrate the effec- 
tiveness of the algorithms and analytical tools contained therein. 
The model has a format similar to the Cucker-Smale ones. The motion of each 
bird's velocity described as follows: 

E J' 1^1 + E {\\^,-x.r-d){x,~x,) 

3 V formation 

For purposes of this paper, we will consider flve birds flying in a V-like forma- 
tion in the context of three different models that have been put forth by various 
authors. We will consider the motion of the birds over time and determine if 
equilibrium will be attained, and if that will lead to long-term stability, based 
on these models. 

We will assume that the birds fly in a coordinate plane parallel to the ground. 
Further, we will consider flve birds cci,..., 2:5, and we will deflne their initial po- 
sitions as follows: 

Xi = (-2,2), 2:2 = (-1,1), a;3 = (0,0), X4 = (-l,-l) 2-5 = (-2, -2). 

Further, we will assume that the initial velocity of the birds is translation to 
the right at velocity (1, 0) so they satisfy: 

= + (t, 0) where j = l,...,5 
We will call this the V evolution. 
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3.1 CFP model 



As stated above, this model considers a force of attraction and a force of re- 
pulsion. Each animal is both attracted and repulsed by the other birds by 
predetermined parameters. 

We want the V evolution to solve the equations as described below: 

\ " {Xj — Xi) 



j(^) Fd ^ ^ (^Xj Xi^ Fr / ^ 



■x^f 



In these equations, Fa and Fr are two positive constants. As we will be setting 
Xiif) to zero, we will need to determine Fa and Fr explicity. 
The region in the last equation depends on parameters a_R (angle of anisotropic 
repulsion regions), Rsr (radius of circular repulsion region) and a a (angle of 
anisotropic attraction region). 

We will consider two cases: 

case \ ar — 90, Rsr — and aa = 179. 

case 2 ar = 360, i?^,. and aa ~ 360. 

Let us first address case 1. Then, animal x^ will not consider the other birds, 
and will therefore keep constant velocity of (1, 0). Let us consider X4. This 
animal will feel both attraction and repulsion towards x^ , and not for the other 
birds. 

Therefore, the equation for the motion of x^ will be defined as follows: 



x^it) = F^{x3 - Xi) - Fy 



{X3 - X4) 
\\X3 - XiW^ 



For Xi above to be a solution to the model, we need the right-hand side of the 
equation to be zero. Taking into account that (xa, X4) = (1, -1) this means 



^^a(l,l) 



Fr 



(i,i) = o. 



thus we need 2Fa=Fr. 

The evolution for each bird in the V formation can be described as: 



(pi(xi,X2, ...jXs). 



where we define 





a;i,i 




Xl 








X2 




X2,l 


X = 


X?, 




X2,2 




X4 








_ ^5 




X5,l 




. 3^5,2 
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The evolution of the whole formation can be written as 
where 













F2 


^3 




F3 








. V5 _ 




Fio 



We will linearize our system around our solution, the V solution. We will 
compute the Jacobian matrix of the vector F of the right hand side of the 
equation above, for i = 1, . . . , 5 along the V solution. Thus F has a total of ten 
components (two for each bird) with each clement of the matrix determined by 
the partial derivatives of F with respect to each component of X as: 

dF, 

For instance, the two rows of the matrix that define the motion of X4 are Fj 

and Fg. 

It follows that 

dFr ^ p _p (3^32 - 3:^42)^ - (3:31 - Xiif) 
dx3,i ° 11 a;3 - 2^4 11"* 

This leads to the value of which is inserted in the matrix. We then compute 
the other values for F7 and Fg using partial derivatives. We perform similar 
calculations for birds xi, X2, and X5, assuming that xi considers birds X2 and 
X3, X2 only considers X3 and X5 considers ^3 and X4. This leads to the following 
matrix: 



-1 


-0.5 


0.5 


0.5 


0.5 

















-0.5 


-1 


0.5 


0.5 





0.5 




















-0.5 


0.5 


0.5 


-0.5 




















0.5 


-0.5 


-0.5 


0.5 






















































































0.5 


0.5 


-0.5 


-0.5 




















0.5 


0.5 


0.5 


-0.5 




















0.5 





0.5 


-0.5 


-1 


0.5 

















0.5 


-0.5 


0.5 


0.5 


-1 



This matrix has eigenvalues of 0, -0.5, -1, -1.5. We can thus write: 

Proposition 1 Consider the CFP model with a,. = 90, = \/2 and aa = 
179. Then the V formation is admissible and not exponentially unstable. 
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Let us now address case 2. In this case, X3 feels both attraction and repulsion 
for X2 and X3 (since we are now considering a circular repulsion region). This 
leads to the following equation: 



However, the model fails to admit V formation as solution. When attempting 
to determine the coefficients for Fa and Fr, it becomes impossible to determine 
constants that will allow the right side of the equation above to be set equal to 
zero. 



3.2 The Cucker-Smale model 

For this model, we will apply similar computations and construct a similar 
matrix to the last model, this matrix will now be 20 by 20 instead of 10 by 10, 
as each bird will have four components in the matrix, two for .T,;(t) and two for 

V^it). 

Since Xi (t) = Vi for each bird, the rows of each matrix that correspond with Xi 
will take the following form in the matrix: 

"0010" 
1 

when computing the partial derivatives with the four columns xa, Xi2, wa, and 
Vi2. The rows that correspond to Vi(t) will be computed with partial derivatives 
as before. The form for each Vi{t) will look as follows: 

for the bird(s) appropriately chosen. 

We will assume that each bird interacts with the other birds as in Case 
1 for the CFP model. That is, bird 3 is independent of the other birds and 
therefore produces zeroes in the matrix for the four rows corresponding to bird 
3. We compute partial derivatives as in the first model, and this generates the 
following matrix: 
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1 . 



This matrix has eigenvalues of 0, 1, and -1. This suggests that the birds wiU 
not attain stabihty or equiUbrium in a simulation of the model. 



Proposition 2 Consider the Cucker-Smale model. Then the V formation is 
admissible and is exponentially unstable. 



4 The Olfati-Saber model 

This model is quite similar to the Cucker-Smale ones, thus we will only indicate 
the main differences. The sum of terms in the derivative of velocity, is only taken 
for birds interacting, so we will consider the birds interacting as in Cucker-Smale 
model {xi only considers X2, X2 only considers x^, x^ is independent, etc.) 

We will define d ~ \pl so that the birds will interact as described before. 
Wc compute partial derivatives as before and compute the following matrix: 

00 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

20 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

00 0000 000 

00 100 000 

00 0010 000 

04-\/2 -2 20 00 

-2-2 4-^20 2 00 

00 0000 010 

00 0000 001 

V2 - 4 2 -2 4 -\/2 -2 2 

2 V2-4 -2 -2 4-V2 2 
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This matrix generates eigenvalues of -1, -0.356, -1.644, -1.363, 3.363, 2.259 
and 0. 

Proposition 3 Consider the Olfati-Saber model. Then the V formation is ad- 
missible and is exponentially unstable. 



5 Conclusions 

We considered three recent models proposed to describe the behaviour of multi- 
agent systems, in particular animal groups. A simple V like formation was 
tested on each model in terms of admissibility of evolution and stability. The 
main result is that V formation are exponentially unstable for two models. This 
indicated the need of further modeling efforts in order to appropriately address 
the evolution of birds flocks. 
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1 Introduction 

Group behavior has received much attention as a test case of self-organization. 
There has been much written in recent years to investigate interactions within 
groups of agents. These agents can be animals moving in an interactive way, such 
as birds, but can also refer to situations such as people driving in traffic. The 
models that describe these interactions are able to reproduce different structures 
and patterns relating to the movement and interaction of the agents involved. 
The advantages and necessities of this type of analysis in any complex biolog- 
ical, technological, economic, or social system are important and far-reaching. 
Each model that we will discuss describes interaction between agents. 
Focusing on animal groups, such systems received attention not only from biol- 
ogists and ethologists, but also from physicists, mathematicians, and engineers. 
This interest has produced a huge amount of literature, which is well docu- 
mented and reviewed: see, for instance, [7, 8, 11] and [3]. Loosely speaking, the 
basic idea behind these works is that complex collective behavior arises from 
simple interactions among close animals. 

This paper considers a simple case of five birds in V formation, flying at 
constant velocity. The aim is to study the feasibility of this solution and its 
stability for three models proposed in the literature. 

The first model we will analyze was proposed in [4]. It is an agent-based model, in 
the sense that each agent is singularly considered. The model is also leaderless, 
meaning that all agents behave following the same set of rules and their behavior 
is not imposed by others. The model is purely based on attraction-repulsion 
infractions between group mates. Attraction allows the group to be formed 
and stay tight, while repulsion allows for collision avoidance between group 
mates and keeps them well spaced. This model does not combine attraction 
and repulsion with velocity alignment, unlike many other models that have 
been proposed. The model also assumes that each agent interacts with a limited 
number of group mates that are in a suitable pre-defined sensitivity zone. The 
authors were able to produce cluster formations, line formations and V-like 
formations by altering the attraction and repulsion parameters of their model. 
The second model was proposed by Cucker and Smale in [6] and additional 
results are available in [10]. The model is of second order, in the sense that 
the dynamics is given for position and velocity of each animal. More precisely. 
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each animal adjusts its velocity to the mating ones, which are in its "interaction 
network" . The adjustment is done towards an average value of the various 
different velocities. Therefore this approach is quite related to the theory of 
consensus. 

The third model is again a consensus-type model proposed by Olfati-Saber in 
[9]. In this case, beside the adjustment of velocity towards an average of mating 
birds, there is a term which forces animal to keep a desired distance d. The 
latter is thought as a suitable distance which taken into account both attraction 
among birds and repulsion for avoidance of conflicting trajectories. 

For each model, wc first show that the V formation is an admissible solution 
(for the first model we need to accurately choose parameters). Then we study 
stability of the linearization along the given V formation trajectory. The result 
is that the linearization of the first model gives rise to negative or zero eigenval- 
ues, while the other models presents also positive eigenvalues. We can conclude 
that the second and third model are exponentially unstable. 
In other words, despite the fact that these models can represent V like forma- 
tions, they are not able to represent stable V like formations. This indicates the 
need of further modelling efforts. 

2 Animal groups models 

Let us start with the model proposed in [4] , which we refer to as the CFP model. 
The animals in the model are represented by point particles, which have simple 
continuous-time dynamics. Given N G N, for alH = 1, . . . , and t G M>o, let 
Xi (t) € represent the position of the i-th animal, whose evolution is described 
by the differential equation 



where x{t) is the vector {xi{t), . . . ,XN{t))- 

The velocity Vi{x) is the sum of two contributions, expressing the effects of 
attraction and repulsion. 



In more detail, each of these contributions depends on the relative position of 
the other animals. 



Xi{t) 



(1) 



v,{x) = vt{x) 






The following definitions have been used. 
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• The function /a : M>o M>o (resp., /r : R>o — > K>o) describes how each 
animal is attracted (resp., repelled) by a neighbor at a given distance, 
assuming || • || denotes the Euclidean norm in R^. 

• The attraction neighborhood (resp., the repulsion neighborhood TZ^) 
is the set of the n animals closest to the i-th one, which are inside the 
attraction (resp., repulsion) sensitivity zone. 

The above model is very general, and we need to specialize it by choosing 
the interaction functions and the shape of the sensitivity zones. We make the 
following assumptions. 

Al. The functions /a and /r are assumed to be 

/adkj - X.,\\) = F^\\Xj - X.,\\ , M\\Xj - X^\\) = ^ r-, 

\\Xj Xi\\ 

where i^a and are two positive constants. 

A2. The sensitivity zones are depicted in Figure 1 and illustrated as follows. 
Let the center point be the animal's position, and let the horizontal axis 
(arrow-headed) represent the direction of motion. Attraction is active in 
a frontal cone whose width is given by the angle aa & (0, 360°] (dashed 
line). Repulsion is active both inside a disk of radius i?sr > {short-range 
repulsion) and in a frontal cone of width € (0°, 360°] (sohd hue). We 
stress that the sensitivity zones do not necessarily coincide with the visual 
field of the animal. They rather represent the zones which attraction and 
repulsion are focused on. 

A3. The speed of each animal \\vi\\ is bounded from above by a constant Wmax- 

The above assumptions result in the system 

i.it) -F^Y. {X, - xO - 1^^-^- (2) 

Some remarks are in order. 

Rl. The definition of the interaction neighborhoods and TZ^ allows to have 
an a priori bound on the number of effective neighbors, and therefore 
on the sensing and "computational" effort which is required for each an- 
imal. This fact, which copes with animals' intrinsic limitations, has been 
experimentally observed in biology, for fish [1] and birds [2]. The latter 
paper calls this neighborhood definition topological, as opposed to metric 
definitions, based on distance only. 

R2. Our assumption of unbounded sensitivity regions does not intend to imply 
that animals sensing capabilities extend on an unlimited range, but rather 
that group dynamics happen in a relatively small area. 
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Figure 1: The shape of the sensitivity zones. 



R3. If Q!a = ar = 360° the parameter Rgr has no effect, and the interaction is 
completely isotropic as in the simulations presented in [2, 5]. Note that, 
even if there is no preference for any specific direction, the limitation of 
the number of considered neighbors makes the interactions not reciprocal, 
i.e. the fact that the z-th animal interacts with the j-th does not imply 
that the j-th interacts with the i-th. 

R4. By specializing the functions /a and /r, one can obtain various interaction 
models. For instance, it is natural to assume that /r be decreasing and 
lims_j.oo /r(s) = 0. Similar considerations are valid also for some features 
introduced in other models, such as a neutral zone around animals or 
a hierarchical decision tree which allows the repulsion force to have the 
priority over the attraction force. 

The second model we will analyze is the Cucker-Smale model. It is inspired 
by the Cucker-Smale flocking model, and attempts to analyze animal behavior 
under hierarchical leadership. The goal of this flocking study is to be able 
to interpret, model, analyze, predict, and simulate various flocking or multi- 
agent behavior. The Cucker-Smale flocking model includes an interaction weight 
function which will allow us to analyze the birds in our formation by defining 
the weight function appropriately. In a model of hierarchical leadership, the 
agents are ranked in some way and are therefore not independently motivated 
in their movement - that is, lower-rank agents are led by agents of higher ranks. 
Using these assumptions, the paper goes on to discuss various flocks using these 
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assumptions. 

The Cucker-Smale flocking model which is given by the autonomous dynamical 
system: 

\ Vi{t) = EjeL(j) - ""0: i = 1, . . . t > 

This model differs from the previous model in that the motion in each bird is 
defined in terms of Xi{t) and Vi{t). In other words, this is a second order model, 
assuming acceleration driven by "social forces". In this case L[i) denotes the 
subgroup of animals that directly influence agent i. 

The third model we will analyze is from [9] . It presents a theoretical frame- 
work for design and analysis of distributed flocking algorithms. In this paper, 
three flocking algorithms are presented: two are for free-flocking and one for 
constrained flocking. As in the Cucker-Smale model, this model compensates 
for velocity alignment. Their first algorithm generally leads to fragmentation, 
whereas the second and third algorithms both lead to flocking behavior. How- 
ever, unlike the second model, this paper asserts that flocks do not need leaders. 
The paper also includes various simulation results that demonstrate the effec- 
tiveness of the algorithms and analytical tools contained therein. 
The model has a format similar to the Cucker-Smale ones. The motion of each 
bird's velocity described as follows: 

3 V formation 

For purposes of this paper, we will consider flve birds flying in a V-like forma- 
tion in the context of three different models that have been put forth by various 
authors. We will consider the motion of the birds over time and determine if 
equilibrium will be attained, and if that will lead to long-term stability, based 
on these models. 

We will assume that the birds fly in a coordinate plane parallel to the ground. 
Further, we will consider five birds cci,..., 2:5, and we will define their initial po- 
sitions as follows: 

11 = (-2,2), X2 = (-l,l), 2-3 = (0,0), X4 = (-l,-l) 2-5 = (-2, -2). 

Further, we will assume that the initial velocity of the birds is translation to 
the right at velocity (1, 0) so they satisfy: 

Xi(t) — Xi^Q + {t,Q) where i = l,...,5 
We will call this the V evolution. 
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3.1 CFP model 



As stated above, this model considers a force of attraction and a force of re- 
pulsion. Each animal is both attracted and repulsed by the other birds by 
predetermined parameters. 

We want the V evolution to solve the equations as described below: 

{Xj — Xj) 



jell" 



\\Xj-Xi\\^' 



In these equations, Fa and are two positive constants. As we will be setting 
Xi{t) to zero, we will need to determine Fa and Fr explicity. 
The region in the last equation depends on parameters (angle of anisotropic 
repulsion regions), R^r (radius of circular repulsion region) and a a (angle of 
anisotropic attraction region). 



We will consider two cases: 

case 1 Ur — 90, Rsr = and = 179. 

case 2 Qfr = 360, R^r ~ and aa ~ 360. 

Let us first address case 1. Then, animal 3:3 will not consider the other birds, 
and will therefore keep constant velocity of (1, 0). Let us consider X4. This 
animal will feel both attraction and repulsion towards x^ , and not for the other 
birds. 

Therefore, the equation for the motion of x^ will be defined as follows: 



x^it) = F^{x3 - X4) - F,. 



{X3 - X4) 
\\X3 - XiW^ 



For X4 above to be a solution to the model, we need the right-hand side of the 
equation to be zero. Taking into account that (0:3, X4) = (1, -1) this means 

F„(i,i)-|:(i,i) = o, 

thus we need 2Fa=Fr. 

The evolution for each bird in the V formation can be described as: 



(pi{xi,X2, —,X5). 



where we define 









Xl 




Xl,2 




X2 




X2,l 


X = 


X3 




X2,2 




X4 












X5,l 
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The evolution of the whole formation can be written as 
where 



<y5l 










F2 


^3 




F3 








. '^5 




Fio 



We will linearize our system around our solution, the V solution. We will 
compute the Jacobian matrix of the vector F of the right hand side of the 
equation above, for i = 1, . . . , 5 along the V solution. Thus F has a total of ten 
components (two for each bird) with each element of the matrix determined by 
the partial derivatives of F with respect to each component of X as: 

dF, 

For instance, the two rows of the matrix that define the motion of X4 are Fj 

and Fg. 

It follows that 

dFr ^ p _ p (3:32 - 3:42)^ - (3:31 - 3:41)'^) 
dxz^i " II xz - Xi |j4 

This leads to the value of ^, which is inserted in the matrix. We then compute 
the other values for F7 and Fg using partial derivatives. We perform similar 
calculations for birds xi, X2, and X5, assuming that xi considers birds X2 and 
X3, X2 only considers X3 and 3:5 considers 3:3 and 3:4. This leads to the following 
matrix: 



-1 


-0.5 


0.5 


0.5 


0.5 

















-0.5 


-1 


0.5 


0.5 





0.5 




















-0.5 


0.5 


0.5 


-0.5 




















0.5 


-0.5 


-0.5 


0.5 






















































































0.5 


0.5 


-0.5 


-0.5 




















0.5 


0.5 


0.5 


-0.5 




















0.5 





0.5 


-0.5 


-1 


0.5 

















0.5 


-0.5 


0.5 


0.5 


-1 



This matrix has eigenvalues of 0, -0.5, -1, -1.5. We can thus write: 

Proposition 1 Consider the CFP model with a,. = 90, Rs,. = and aa = 
179. Then the V formation is admissible and not exponentially unstable. 
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Let us now address case 2. In this case, X3 feels both attraction and repulsion 
for X2 and X3 (since we are now considering a circular repulsion region). This 
leads to the following equation: 



X3{t) = Fa{{x2 - X3) + {X4 - X3)) - Fr{ 



{X2 - X3) 



+ 



{X4 - X3) 



)■ 



\\X2-X3r 



\\X4 - X3I 



However, the model fails to admit V formation as solution. When attempting 
to determine the coefficients for Fa and Fr, it becomes impossible to determine 
constants that will allow the right side of the equation above to be set equal to 
zero. 

3.2 The Cucker-Smale model 

For this model, we will apply similar computations and construct a similar 
matrix to the last model, this matrix will now be 20 by 20 instead of 10 by 10, 
as each bird will have four components in the matrix, two for a;,;(t) and two for 



Since Xi{t) = Vi for each bird, the rows of each matrix that correspond with Xi 
will take the following form in the matrix: 



when computing the partial derivatives with the four columns xn, Xi2, ua, and 
Vi2. The rows that correspond to Vi(t) will be computed with partial derivatives 
as before. The form for each Vi{t) will look as follows: 



for the bird(s) appropriately chosen. 

We will assume that each bird interacts with the other birds as in Case 
1 for the CFP model. That is, bird 3 is independent of the other birds and 
therefore produces zeroes in the matrix for the four rows corresponding to bird 
3. We compute partial derivatives as in the first model, and this generates the 
following matrix: 



V^{t). 



10 
1 



8 



" 





1 


















































" 











1 
























































— 1 











1 


















































— 1 











1 
























































1 






























































1 
























































— 1 











1 


















































— 1 











1 
























































1 






























































1 




























































































































































































1 






























































1 












































-1 











1 


















































-1 











1 




































































1 






























































1 












































-1 











1 





. 












































-1 











1 . 



This matrix has eigenvalues of 0, 1, and -1. This suggests that the birds wih 
not attain stabihty or equiUbrium in a simulation of the model. 



Proposition 2 Consider the Cucker-Smale model. Then the V formation is 
admissible and is exponentially unstable. 



4 The Olfati-Saber model 

This model is quite similar to the Cucker-Smale ones, thus we will only indicate 
the main differences. The sum of terms in the derivative of velocity, is only taken 
for birds interacting, so we will consider the birds interacting as in Cucker-Smale 
model {xi only considers X2, X2 only considers X3, X3 is independent, etc.) 

We will define d = V2 so that the birds will interact as described before. 
Wc compute partial derivatives as before and compute the following matrix: 









1 


















































1 






































V2 - 4 


-2 





4 - V2 


-2 


2 
































2 


V2 - 4 





-2 


-2 


4 - \/2 





2 












































1 


















































1 






































^2 -4 


2 


-2 





4-^2 


-2 


2 





























2 


\/2 - 4 





-2 


-2 


4 - \/2 





2 
























































































































































































































































1 


















































1 


























V2-4 


2 


-2 





4 - \/2 


-2 


2 





























2 


V2 - 4 





-2 


~2 


4 - V2 





2 








































































































































2 


-2 









































2 


^2 - 4 





-2 
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4 - \/2 -2 2 

-2 4 - ^2 2 
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This matrix generates eigenvalues of -1, -0.356, -1.644, -1.363, 3.363, 2.259 
and 0. 

Proposition 3 Consider the Olfati-Saber model. Then the V formation is ad- 
missible and is exponentially unstable. 



5 Conclusions 

We considered three recent models proposed to describe the behaviour of multi- 
agent systems, in particular animal groups. A simple V like formation was 
tested on each model in terms of admissibility of evolution and stability. The 
main result is that V formation are exponentially unstable for two models. This 
indicated the need of further modeling efforts in order to appropriately address 
the evolution of birds flocks. 
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